Quenched lattice calculation of the vector channel B — > D*l.i> decay rate 
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We calculate, in the continuum limit of quenched lattice QCD, the form factor that enters the 
decay rate of the semileptonic decay B — > D*tv. By using the step scaling method (SSM), previously 
introduced to handle two scale problems in lattice QCD, and by adopting flavor twisted boundary 
conditions we extract F B ^ D (w) at finite momentum transfer (w > 1) and at the physical values 
of the heavy quark masses. Our results can be used in order to extract the CKM matrix element 
Vcb by the experimental decay rate without model dependent extrapolations. The value of V c b 
agrees with the one obtained from the B —* Diu channel and makes us confident that the quenched 
approximation well applies to these transitions. 



I. INTRODUCTION 



Physics beyond the Standard Model may show up in 
the hadronic flavor sector of the theory and could be re- 
vealed by measuring independently the different entries 
of the Cabibbo-Kobayashi-Maskawa [l|, 0] matrix and by 
looking for deviations from unitarity (unitarity triangle 
analysis, UTA) [3|. On the theoretical side, a percent 
relative accuracy on the hadronic matrix elements enter- 
ing the different variants of the UTA's 4] is needed and 
lattice QCD calculations may eventually provide the re- 
quired non-perturbative precision. The study of semilep- 
tonic decays of heavy-light mesons mediated by flavor 
changing AF = 1 hadronic currents gives direct access 
to CKM matrix elements and is particularly convenient 
from the point of view of lattice QCD calculations. The 
form factors entering the decay rates are dimensionless 
quantities and do not inherit uncertainties from the scale 
setting procedure. Moreover they can be expressed in 
terms of matrix elements undergoing finite and multi- 
plicative renormalization and symmetry arguments can 
be used to constrain numerical results (Ademollo-Gatto 
theorem, Luke's theorem, spin-symmetry in the heavy 
quark limit). Finally, matrix elements involve a single 
hadron both in initial and final states thus avoiding com- 
plications due to final state interactions. 

In two previous papers (fj, [6j we have calculated, in 
the quenched approximation of QCD, the form factors 
entering the decay rate of the process B — > Dlv at non 
vanishing momentum transfer and, comparing with ex- 
periment, we have extracted the CKM matrix element 
V c b- In this work we extend our quenched study of 
heavy-light mesons semileptonic decays by calculating 
the form factors entering the decay rate for the pro- 
cess B — > D*tv. At vanishing momentum transfer the 
study of pseudoscalar to vector semileptonic transitions 
requires the calculation of a single form factor with re- 
spect to the two required in the case of pseudoscalar 
to pseudoscalar transitions. At finite momentum trans- 
fer the study of pseudoscalar-vector transitions requires 
the calculation of a particular linear combination of four 



different form factors, F B ^ D (w). Here we compute 
pB^D ^ w ^ j n j.jjg m t erva j i < w — VB . VDt < i_i 

thus allowing the extraction of V c b without extrapolat- 
ing experimental data to zero recoil. Although in this 
case the extrapolation is much less dramatic than in the 
pseudoscalar-pseudoscalar channel an important source 
of systematics can be avoided by using our results. 

As in refs. [E|, [y] we make use of the Step Scaling 
Method [7| devised to reconcile large quark masses with 
adequate lattice resolution and large physical volumes 
and successfully applied also to the determination of 
heavy quark masses and decay constants @, M, Ell • We 
obtain results at non zero momentum transfer with good 
accuracy by enforcing special boundary conditions on the 
quark fields, called flavor twisted [ll[, that shift by an 
arbitrary amount the discretized set of lattice momenta 
(see also CHllEi). 

Our results are not the final ones since they have been 
obtained within the quenched approximation. Quench- 
ing introduces a systematic error that it is hard to quan- 
tify (if not impossible) but allows us, in view of a fu- 
ture unqucnchcd calculation, to check all the remaining 
systematics (heavy quark methodology, continuum and 
chiral limits, etc.) and to discuss some technical issues 
related to the choice of interpolation operators for vector 
mesons carrying non vanishing spatial momenta to be 
used within the Schrodinger Functional formulation of 
lattice QCD. Furthermore precise results for F B ^ D (w) 
at w > 1 are presently missing even in the quenched 
approximation. 

To estimate the validity of the quenched approxima- 
tion we calculate the ratio of physical decay rates between 
vector and pseudoscalar final states. This ratio is indeed 
a purely QCD observable, independent from the value 
of V c b, and can be directly compared with experiment. 
Our results agree with the measured ratio within exper- 
imental errors thus indicating that residual unquenched 
corrections are likely within the current experimental un- 
certainties. 



II. FORM FACTORS AND DECAY RATE 

The semileptonic decay of a pseudoscalar meson into a 
vector meson is mediated by the weak V — A current. The 
relevant matrix elements can be parametrized in terms of 
four form factors. Among possible parameterizations we 
choose the following one 
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where we have used the greek letters n,v,p,a to indicate 
Lorentz indices, Mpy are the masses of the pseudoscalar 
(Mp) and vector (My) mesons, vpy = Pp,v/Mp t v their 
4-velocities, e^ pa is the totally antisymmetric tensor in 
four dimensions (e 0123 = 1) while e p is the polarization 
vector of My, 
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(2) 



Q = l 



The form factors depend upon the masses of the initial 
and final particles and upon w = vy ■ vp 
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1 < w < (Mf, + M v )/2MpM v 



(3) 



In the case where Mp is the B meson mass and My is 
the D* meson mass the maximum value of w is around 
1.5. 

The differential decay rate of the process B — ► D*£v, 
in the case of massless leptons, is given by 



dT 



B^D*lv 



dw 



= IK h | 2 ^Mi(l-r)Vx 



48tt 3 ' 



a/io 2 - 1(1 + wf\(w) F 

where we have defined r — My /Mp and 

1 — 2wr + r 2 



? B^D' 



(w) 



(4) 



*» 



(1-r)* 



A(w = 1 + t 2 (w) 

w + 1 



(5) 



m, 



pfl ->D*, v , , , H$(w)+Hl(w)+Hl(w) 
F (w) = h Al (w)\ 



AH 



(6) 



with 
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is the non perturbative input needed to extract V c b by 
the measurement of the decay rate. 

By noting that at zero recoil F B ~* D (1) is identically 
equal to h Al (l) the complexity of the theoretical calcu- 
lation can be considerably reduced since, in this partic- 
ular case, a single matrix element with initial and final 
particles both at rest is needed instead of the four ma- 
trix elements at non vanishing momentum transfer re- 
quired to solve the full system of eqs. |T]) with respect 
to hy,Ai ■ For these reasons previous lattice studies have 
been devoted to the calculation of F B ^ D (w) at zero re- 
coil only, where it can be extracted with good statistical 
accuracy both in the quenched approximation [lq and 
in the rif — 2 + 1 unquenched theory [15] . On the other 
hand, it is not possible to measure directly the decay rate 
at zero recoil because of the presence of the kinematical 
factor (w — l) 1 / 2 in eq. Q) and experimental data at 
w = 1 are obtained by extrapolation. The systematics 
introduced by this extrapolation is much less dramatic 
with respect to the case of the decay B — ► Dlv where 
the kinematical suppression goes like (w — l) 3 / 2 but can 
be nevertheless eliminated. In the following we calculate 
pB^D j-^-j j n j.jjg ran g e 1 < ui < 1-1 that includes values 
of w where experimental data are directly available with 
good precision. 



III. SCHRODINGER FUNCTIONAL 
CORRELATORS 



We have carried out the calculation within the Oia) 
improved Schrodinger Functional formalism [171 IX 81 ] 



with vanishing background fields. The choice of the 
Schrodinger Functional regularization is particularly con- 
venient to perform simulations on small physical volumes 
(see section HVl) because Dirichelet boundary conditions 
in the time direction provide an infrared regulator that 
allows the simulation of massless quarks. At the same 
time, the extraction of physical matrix elements involv- 
ing vector mesons at non vanishing spatial momenta from 
Schrodinger Functional correlators requires some addi- 
tional care with respect to the case of quark fields satis- 
fying periodic boundary conditions in the time direction. 



Indeed the two boundary projectors act on 7^ by killing 
its time component and by introducing a mixing with the 
tensor a ^ = i[j°,j' 1 ]/2, 
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p = 
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i(7 l -kr w ), n = i = 1,2,3 



(12) 



which, in compact notation, is equivalent to the gamma 
matrix combination appearing in eq. (jlip 

2V - 2tV'° + h° VI 
P+^P- = =-? Z7 9 +17 ,7 ] (13) 



A. Boundary and bulk operators 



In defining interpolation operators of meson states we 
need to distinguish between bulk fields tj}{x), boundary 
fields ((x) living at x — 

p+C(x) = i±™ vM = o 



and boundary fields C( x ) living a t xq = T 
P-C'W = l —^^{T,x)=0 



(8) 



The matrix elements of V M between the vacuum and a 
vector meson state, entering the spectral decomposition 
of two and three point correlation functions, can be thus 
parametrized as 

(Q\V»\M V ) = p v [e£ - e° a g°»] + PT [e»v v - e° a v v ] 

(14) 

The projection of the wall sources on the physical states 
entail unknown non-perturbative wave functions, pv and 
PT-, that cancel out exactly in our choice of ratios of corre- 
lation functions discussed in the next subsection. Pseu- 
doscalar mesons do not carry polarization indexes and 
the effect of the boundary conditions can be reabsorbed 
into a redefinition of the wave function, 



C'(aQP+ - ^(T,f)i±^ = (9) 

Different quark flavors will be distinguished, if needed, 
by using explicit indexes. External momenta have been 
set by using flavor twisted b.c. for the heavy flavors. In 
particular we have used 
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with different values of 9 for the different heavy quarks 
and ordinary periodic b.c. in the other spatial directions 
and for the light quarks. 

A generic meson state on the boundaries can be ex- 
pressed as a bilinear field operator acting on the vacuum. 
In the case of vector mesons, the Schrodinger Functional 
boundary conditions select a particular combination be- 
tween the two possible choices V^{x) = , ip{x') r ) ,1 rp{x) and 
T 0/i (x) = ^(x)7 O 7 /J '0(x). On the wall at xq = we have 
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(15) 



The relations above hold at xq — but the same argu- 
ments can be repeated for the boundary operators on the 
wall at xq = T that in the following will be called P' and 
V"\ 

For later use, we also define the improved bulk opera- 
tors 

A»(x Q ) = A»(x) + ac A ^-^P{x) 



V(x) = V(x) + ac v dv \ d * v T^(x) (16) 



a" 



E ^(°. y)- — \ > (°' z ) 



(11) 



where P(x) = ip(x) r YBip(x) is the pseudoscalar density, 
A ll {x) = iplx^^j^iplx) the axial current and the im- 
provement coefficients ca and cy have been taken from 
refs. [UllliH. 



B. Quark masses 



Particle states are fixed by timing the values of the 
corresponding quark masses. We use three definitions of 
renormalization group invariant quark masses differing at 
finite lattice spacing by terms of order 0(a 2 ). We have 
calculated the following two point correlation functions 



f r A r (x ) = J^(P rr A° rr (x)) 
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/£(*„) = -J2(PrrPrr(x)) 



(17) 



and defined 
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(18) 
(19) 



where a is the lattice spacing, k r is the hopping parameter 
of the r quark and k c is the critical value of the hopping 
parameter. 

A first definition of renormalization group invariant 
(RGI) quark masses has been obtained by the following 
relation 



= Z M [l + (b A - b P ) am b r ] 



AWI 



(20) 



The combination b A — bp of the improvement coefficients 
of the axial current and pseudoscalar density has been 
computed non-perturbatively in [22, [2J] . The factor Zm 
is known with very high precision in a range of inverse 
bare couplings that does not cover all the values of j3 used 
in our simulations. We have used the results reported in 
table 6 of ref. [23| to parametrize Zm in the enlarged 
range of (3 values [5.9,7.6]. A second definition of RGI 
quark mass has been obtained by the relation 



= Z M Z [l + b m am b ] 



amJ\ m_ 



(21) 



where the improvement coefficient b m and the renormal- 
ization constant Z has been also taken from refs. [22|, |2J] . 
An additional definition of RGI quark masses has been 
obtained by using improved lattice derivatives in eq. (fT5| . 



C. Three point correlators 

The matrix elements in eqs. (JXJ) can be calculated on 
the lattice by building suitable ratios of the following 



three point correlation functions 

(PVP)? f (x ,p Pl ,f Pf ) = ^S^W^') 
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(22) 

In the previous equations we have explicitly indicated 
heavy (z, /) and light (/) flavor indices and we have im- 
plicitly defined 
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(23) 



Let us now consider the normalization factors 

Nif( x o,m,pv f ) 

= y/(PVP) % (x , m , p Pi ) (VVV) f* (x , p Vf , Pv f ) 

(24) 

Because of our choice of spatial momenta (having non 
vanishing components along the direction 1 only) the nor- 
malization factors A 1 will in general be different from 
the remaining ones, A 2,3 . By assuming single state dom- 
inance, by relying on the conservation of the vector cur- 
rent and by using the form factors definition, eqs. |T]), 
and the completeness relation of the polarization vectors, 
eq. ^, we have 
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V v p v v\ N 2 (T/2,p P ,p v ) 
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V W 2 — lfly 



(25) 



The extraction of h Al from the first of the previous re- 
lations is straightforward. In view of the calculation of 
the decay rate (see eqs. ((U, dU and ©), we do not re- 
move the factor \Jw 2 — 1 from the second equation above. 
More precisely, we measure on the lattice the following 
combinations 
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three point correlation functions by the following ratios 

l(VVV)U l (T/2,p Vz ,p v ,) 
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that enter directly the definition of F B ^ D (w). In the 
case of the remaining form factors, h A2 and h Aa , we fol- 
low a similar path. By defining 



PP^ = (PVP)l(T/2,p Pt ,p P J 
E Pi (PVP)%(T/2,p Pi ,p Pi ) 



(32) 



This way we have been able to define all the form factors 
entirely in terms of three point correlation functions and 
to keep under control statistical fluctuations also at non 
vanishing momentum transfer. 
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we measure on the lattice the following combinations of 
correlation functions 



X 2 = (w- 1) 



h A2 v v C'-v v C° 

hZ- — c^ l (29) 



X 3 = (w-l)-^ = 
h Al 



Aa v p C° ~ VpC 1 
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that provide the combinations of form factors entering 
the definition of F B ^ D (w). 

A key point in reducing statistical fluctuations in the 
extraction of the different combinations of form factors is 
the definition of the 4- velocities, and consequently of w — 
vp ■ w, entering eqs. (|2"6")l . (f2U)) and (|3"D|) . A particularly 
convenient choice is the definition of velocities in terms of 



IV. THE STEP SCALING METHOD 



The idea behind the SSM is to split the complex- 
ity of the calculation of quantities depending upon two 
largely separated energy scales into several calculations 
performed on different physical volumes. The small vol- 
ume calculations are needed in order to resolve the dy- 
namics of the heavy quarks without recurring to any ap- 
proximation but introducing, at intermediate stages, fi- 
nite volume effects (FVE) 

O(physioal) = 0{finite volume) x FVE (33) 

The finite volume effects are subsequently accounted for 
by performing simulations on progressively larger vol- 
umes. The success of this strategy depends on the details 
of the problem and hence on the possibility of computing 
the finite volume observable, finite volume effects and 
their product with smaller errors and systematics with 
respect to the ones that would be obtained by a direct 
calculation. The strength of the method is a great free- 
dom in the definition of the observable on finite volumes 
provided that its physical value is recovered at the end 
of the procedure. 

In the calculation of heavy-light meson observables the 
two scales are the mass of the heavy quarks (b,c) and the 
mass of the light quarks (u,d,s). Here we consider the 
form factor F p ^ v (w) as a function of w, the volume, 
L 3 x T, and identify heavy meson states by the corre- 
sponding RGI quark masses that in the infinite volume 
limit lead to the physical meson spectrum [9(. 

First we compute the observable F B ~* D (w; Lq x Lq) 
on a small volume Lq x Lq, T = L, Lq ~ 0.4 fm, which 
is chosen to accommodate the dynamics of the 6-quark. 
The observable on the finite volume is defined as ex- 
plained in the previous section, i.e. by relying on single 



state dominance of the different correlation functions. Of 
course this is not true at xq = L /2 ~ 0.2 fm but, letting 
xq to scale proportionally to the volume in the finite size 
scaling iteration, becomes true by removing finite volume 
effects. 

A first portion of finite volume effects is removed by 
evolving the volume from Lq x Lq to L\ x L\, L\ = 2Lq, 
by the ratio 



_P—>D* 



(w;L ,Li) 



F 



P-*D* 



{w\L\ x L x ) 



Fr->°*{w;LlxL ) 



(34) 



The crucial point is that the step scaling functions 
are calculated by simulating heavy quark masses mp 
smaller than the 6-quark mass. The physical value 
a B ^ D (w; Lq, Li) is obtained by a smooth extrapolation 
in 1/mp that relies on the HQET expectations and upon 
the general idea that finite volume effects, measured by 
the tr's, are almost insensitive to the high energy scale. 

The final result is obtained by further evolving the 
volume from L\ x L\ to L\ x T 2) with L 2 = 4Lo an d 
T 2 = 3L 2 /2 by calculating a second step scaling function 



L A 
L B 
LoC 



L a 
L b 
Loc 

L X A 
Lj_B 
LxG 



L\a 
L ± b 
L x b 

L 2 A 
L 2 B 
L 2 C 



a 



L 3 xT 



N t 



7.6547 
7.4082 
7.2611 



7.0203 
6.7750 
6.4956 

7.0203 
6.7750 
6.4956 



6.4956 
6.2885 
6.0219 

6.4956 
6.2885 
6.0219 



32 3 x 32 
24 3 x 24 
20 3 x 20 
16 3 x 16 
12 3 x 12 
8 3 x 8 

32 3 x 32 
24 3 x 24 
16 3 x 16 



16 3 x 16 
12 3 x 12 

8 3 x 8 

32 3 x 48 
24 3 x 36 
16 3 x 24 



cnjg 



433 
298 
128 



293 
640 
1600 

326 
243 
120 



1051 
3168 
720 

152 
150 
200 



TABLE I: Table of lattice simulations. 
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and by the following identity 



F 



B^D* 



(w;L 3 2 x3L 2 /2) 



F 



B^D* 



[w; LI x L a 



x <j B ^ D \w\L ,L x ) 

x a B ^ D '(w;L u L 2 ) (36) 



LATTICE SIMULATIONS 



Small volume 



Simulations on the small volume Lq x Lq have been 
carried out at three different values of the lattice spac- 
ing corresponding to the entries LoA, LqB and LqC in 
TABLE |I1 The beta values have been deduced from 
refs [25|, [2y, [27| and have been matched by fixing the 
value of the renormalized strong coupling constant in the 
Schrodinger Functional scheme. The physical extension 
of the volume is Lq ~ 0.718ro where tq ~ 0.5 fm 28]. 

On such a small volume we have simulated massless 
light quarks and we do not need to discuss the systemat- 
ics associated to chiral extrapolations. 

In FIG. [T] we show the continuum extrapolation of 
pB^D ^y. ^3 x £ ) f or th e seven different values of w 
that have been simulated. Different fits correspond to the 



three different definitions of quark masses introduced in 
section IIIII the extrapolations must coincide within the 
errors as happens to be in practice and the continuum 
results are defined by a jackknife average. In general 
we combine results of different simulations in big jack- 
knife samples according to the recipe discussed in the 
Appendix A of ref. [29J . 

In FIG. [2] we show the different form factors that en- 
ter into the definition of F B ^ D (w; Lq x Lq) already ex- 
trapolated to the continuum limit. As expected, Xy(w), 
X 2 (w) and X^{w) vanish at zero recoil while hA 1 {w) is 
of order one. The scale of the figure does not allow to 
distinguish the errors that are much smaller than the 
symbols but it allows to appreciate the relative sizes of 
the different form factors. 

Concerning the renormalization factors, our definitions 
of hA 1 {w) and Xy(w) require the knowledge of Za/Zv 
and of (bA ~ by) while, in the case of X 2 (w) and X 3 (w), 
these cancel out in the ratios. Unfortunately the ratio 
Za/Zv has not been computed directly and we have 
used the separate non-perturbative determinations of Za 
and Zy performed in ref. [30j. Also by has been non- 
perturbatively determined in ref. [30] while for bA we 
use its perturbative value (approximation well justified 
at this small values of bare couplings and, a posteriori, 
by looking at the continuum extrapolations of FIG. [1]). 
The overall systematics due to renormalization factors 
is largely accounted for by adding a 0.6% relative error 
to our final result (see ref. [30] ) - Since renormalization 
factors cancel out in the definition of the step scaling 
functions there are no other systematics associated to 
the renormalization procedure in our calculation. 
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FIG. 1: Continuum extrapolations of F B ^ D (w; Lq x Lq) for 
the different simulated values of w. The different fits are 
linear in (a/Lo) and correspond to different definitions of 
RGI quark masses. 



d * 



FIG. 2: Small volume results for ^Ai(w), Xv{w), X2{w) and 
Xs(w). On this scale the errors are much smaller than the 
symbols. 



this step the light quark has been simulated at zero mass. 
The value of the heaviest quark mass simulated on these 
volumes has been halved with respect to the previous 
step in order to have the same order of cutoff effects. 

In FIG. [3] we show the continuum extrapolation of 
fjP^D j-^. £ 0j 2^) for the seven different values of w that 
have been simulated. The pseudoscalar meson state cor- 
respond to a quark having a mass of about half of the 
physical value of the 6-quark mass (mp ~ m\,jT) i.e. the 
heaviest mass simulated on these volumes; similar plots 
could have been shown for the other simulated heavy 
quark masses. Different fits correspond to the three 
different definitions of quark masses introduced in sec- 
tion [nTl the continuum results are defined by the jack- 
knife average of the independent extrapolations. 

In FIG. 0] we can test our hypothesis of the low 
sensitivity of the finite volume effects upon the heavy 
quark mass. The figure shows the step scaling function 
a p ^ D (w; Lq, Li) in the continuum limit and at fixed w 
as a function of the inverse RGI heavy quark mass of the 
pseudoscalar state 1/mp. The RGI heavy quark mass of 
the vector state has been held fixed to its physical value 
m c . The physical step scaling functions are obtained by 
linear extrapolations so mild that the values at rrib differ 
by the simulated ones by a few per millc. 



B. Intermediate volume 



Final volume 



The step toward the intermediate volume L\ x L\, 
L\ = 2Lo, has been performed by calculating the de- 
nominator of eq. (|M|) at three different lattice spacings, 
corresponding to the entries L^a, L$b and Lqc of TA- 
BLE [H and the numerator with the same lattice spacings 
but with twice the number of lattice points per direc- 
tion, entries L ± A, L X B and L X C of TABLE [J Also for 



The last step of the finite size scaling procedure 
has been performed by calculating the denominator of 
eq. (|55|) at three different lattice spacings, corresponding 
to the entries L\a, Lib and L\c of TABLEfl] and the nu- 
merator with the same lattice spacings but with twice the 
number of lattice points along the spatial directions and 
three times the number of points along the time direc- 
tion, entries L 2 A, L 2 B and L 2 C of TABLEQ] The choice 
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FIG. 3: Continuum extrapolations of a p ^ D (w\ Lq, L±) for 
the heaviest quark mass (mp a± rrib/2) simulated on this step 
and for the different simulated values of w. The different fits 
are linear in (a/Lo) 2 and correspond to different definitions 
of RGI quark masses. 



FIG. 4: Continuum step scaling function o p ~^ D (w\Lo,Li) 
as a function of the inverse RGI heavy quark mass of the 
pseudoscalar state, 1/mp, for the different values of w. The 
vertical black lines represent the physical values of the charm 
and bottom quark masses. 



of T = 3L/2 in the numerator of the last step is moti- 
vated by the need of reaching a time extent in physical 
units that justifies the single state dominance hypothesis. 
This can also be compared with the different sequence of 
steps, carried out in our previous papers [3, |g|, reaching 
the same final time extent. 

In FIG. [5] we show the continuum extrapolation of 
a p ^ D (w; Li, L2) for the seven different values of w that 



have been simulated. The pseudoscalar meson state 
correspond to a heavy quark having a mass of about 
a quarter of the physical value of the 6-quark mass 
(mp ~ mfr/4), the heaviest mass simulated on these vol- 
umes; similar plots could have been shown for the other 
simulated heavy quark masses. In this step we have used 
a single definition of RGI quark masses and precisely the 
one of eq. (|20|) . 
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FIG. 5: Continuum extrapolations of a p ~ ,D (w;Li,Li2) for 
the heaviest quark mass (mp ~ m&/4) simulated on this step, 
mi = m 3 and for the different simulated values of w. 



In FIG. Owe show the extrapolations of the continuum 
step scaling function a p ^ D (w; Li, L2) to the physical 
point. The RGI heavy quark mass of the vector state 
has been held fixed to its physical value m c . Also in this 
second step the dependence of the finite volume effects 
upon the heavy quark mass of the pseudoscalar state is 
very mild. The fact that we measure larger finite volume 
effects with respect to the ones obtained on the previous 
step can be explained by assuming that these are mainly 
due to the contribution of excited states to the correla- 
tion functions entering the definition of the lattice matrix 



FIG. 6: Continuum step scaling function a p ^ D (w;Li,Ii2) 
as a function of the inverse RGI heavy quark mass of the 
pseudoscalar state, 1/rap, for the different values of w. The 
vertical black lines represent the physical values of the charm 
and bottom quark masses. Symbols of different colors corre- 
spond to different values of the light quark mass. 



elements rather than to the matrix elements themselves. 
Indeed, in passing from the denominator to the numer- 
ator of eq. (p?5|) the time extent of the final volume is 
enlarged of a factor three, X2 ~ 2.4 fm. The validity 
of this hypothesis is supported by the results of a series 
of simulations of the volume L\ x Li not shown in this 
paper. 
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w 


F B ^ D *(w) 




N f 


reference 


1.000 


0.917(08)(05) 


0.878(10)(04) 





this work 


1.010 


0.913(09)(05) 


0.883(09)(04) 





this work 


1.025 


0.905(10)(05) 


0.891(09)(04) 





this work 


1.050 


0.892(13)(04) 


0.905(10)(04) 





this work 


1.070 


0.880(17)(04) 


0.914(12)(05) 





this work 


1.075 


0.877(18)(04) 


0.916(12)(05) 





this work 


1.100 


0.861(23)(04) 


0.923(16)(05) 





this work 


1.00 


0.913(20)(16) 







[16] 


1.00 


0.924(12)(19) 




2+1 


[15] 



TABLE II: Final results in the continuum and infinite volume 
limits. As a comparison we quote also the results of previous 
lattice calculations of F B ^ D (w) by the Fermilab collabora- 
tion. 



It has not been possible to simulate light quarks at van- 
ishing mass on the final volume and we have computed 
the step scaling function a p ^ D (w;Li,L,2) at three dif- 
ferent values of mi ranging from about 3m s /2 to about 
m s , the physical value of the strange quark mass. Sym- 
bols of different colors in the plots of FIG.[6]show the step 
scaling functions for different values of mi . We do not see 
any appreciable dependence of finite volume effects upon 
mi within the statistical errors and obtain results in the 
chiral limit by a linear fit. The slopes of these fits are at 
least an order of magnitude smaller than their error. 



VI. COMPARISON WITH EXPERIMENT 



FIG. 7: Comparison of \V cb x 10 3 |F S ~* D * (w) obtained in this 
work with experimental data. Lattice data carry only statis- 
tical errors and are normalized at w = 1.075. 




••• 



F(w)/G(w) HFAG fits for PDG 2008 



FIG. 



Comparison of the ratio F B ~' D (w)/G E 



\w) ob- 
tained in this work with the fits to experimental data per- 
formed by the HFAG for the 2008 edition of the PDG. Lattice 
data carry only statistical errors. 



Our final results for F B ^ D (w) are obtained by multi- 
plying the small volume numbers with the two step scal- 
ing functions according to eq. (|36p and are shown in TA- 
BLE |nj The first error is statistical while the second 
comes from the uncertainties associated to the renormal- 
ization factors on the small volume (see discussion in sec- 
tion |TT]) . As a comparison we quote in table TABLE |TT] 
the results of previous lattice calculations obtained by 
the Fermilab lattice collaboration [15|, [jj| at zero recoil. 

FIG. [7] shows the comparison of our lattice data 
with some experimental determinations of \V c t> x 
10 3 \F B ^ D *(w) [1H HI H Hi]: the functional depen- 
dence of experimental data upon w is reproduced by lat- 
tice data within statistical errors also in the quenched 
approximation. The comparison is made by matching 
lattice data with the experimental ones from ref. [3l| at 
w = 1.075 and by obtaining 



\Vcb\ = 3.74(8)(5) x IO- 



CS?) 



where the first error is from theory while the second 
from experiment. This number is in good agreement 
with our previous determination of the CKM matrix el- 
ement, V cb = 3.84(9) (42) x 10 -2 , performed in ref. [1] by 



matching G B ^ D (w), the form factor entering the decay 
rate of the process B — > D£i>, with experimental data 
at w = 1.2. In both the cases V c b has been extracted at 
non vanishing momentum transfer and the corresponding 
uncertainties are not correlated with the errors on the pa- 
rameters of the fits needed to extrapolate measured decay 
rates to zero recoil. 

A more direct comparison between lattice and ex- 
periment can be obtained by introducing the ratio 
F B ^ D * (w)/G B ^ D (w) that does not depend upon V cb . 
In order to apply a jackknife procedure directly to the ra- 
tio of the two form factors, we have computed G B ^ D (w) 
on the same gauge ensambles of F B ^ D (w) by following 
the finite size scaling recursion described in section IIVI 
and by using the definition of G B ^ D (w) discussed in 
refs. p, Q. Our results are given in TABLE HI1 and com- 
pared with experiment in FIG. [5) The black curves have 
been drawn by using the two independent fits 1 of the ex- 
perimental data for F B ^ D (w) and for G B ^ D {w) per- 



1 A direct evaluation from the experimental collaborations of the 
ratio of decay rates would likely lead to a more stringent test. 
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formed by the Heavy Flavour Averaging Group [35| by 
using the parametrizations of ref. [36J]. Our results, red 
points in the figure, fall inside the allowed experimental 
region, i.e. between the two solid black curves. Theo- 
retical errors are smaller than experimental ones that in 
turn are dominated by the uncertainties on G B ^ D (w). 
Also from this parameter-free comparison with experi- 
ments we obtain a strong indication for the validity of 
the quenched approximation for these observables. To- 
gether with the evidence of mild volume effects (if we ex- 
clude those coming from excited states contributions to 
correlation functions) these results may call for a basic 
short distance nature of the form factors. This could jus- 
tify why main unquenching effects are reabsorbed by the 



renormalization procedure of the quenched calculation 
that fixes indirectly the renormalized coupling constant 
from physical quantities. 
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